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ABSTRACT 

Safety analyses of a loss-of-coolant accident in a 
pressurized-water-reactor nuclear power plant require 
modeling the flow through the main-coolant pumps in any one 
of several extreme operating conditions. A semi-empirical 
method of predicting the two-phase performance of centri- 
fugal pumps has been proposed and was applied to normal pump 
operation with good results. In this report the proposed 
method is applied to the reverse-flow, reverse-rotation mode 
of operation, also with good results. 

By making simplifying assumptions, the flow through 
the pump is idealized and expressions for the theoretical 
pump characteristics derived from Euler's equation for the 
mean flow through a turbomachine. General expressions for 
the flow losses in two-phase flow are derived, and the flow 
losses expressed as the difference between the theoretical 
characteristics and the actual performance of the pump. 

The ratio of the losses in two-phase flow to those in single- 
phase flow is shown to be primarily a function of pump 
geometry and fluid void-fraction. Using experimental data 
from single- and two-phase flow tests on a pump of known 
geometry, the ratio of the losses was calculated for various 
conditions of flow and plotted as a function of void fraction. 
The results plotted in this fashion show excellent correlation 
and the function can be used to predict the two-phase 
characteristics of similar pumps. 
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LIST OF SYMBOLS 



A 


flow area 


a 


P 

Ot ~ v 

two-phase flow function = [•= ] — 

i-ce P L 


c 


fluid velocity, absolute 


f tp 


two-phase flow function (eq. 8a) 


g 


gravitational acceleration 


g c 


constant in Newton's Law 


H 


head 


H* 


AH ... - AH . 

head-loss ratio = t- ° ^ AH ° P 

AH ospth ' AH osp 


h 


enthalpy 


K, K 4 , 


, Kg constants defined by eqs. 12, 13, and 15 


• 

m 


mass flow per unit time 


N 


rotor speed, rev/min 


P 


fluid pressure 


Q 


volume flow per unit time 


r 


impeller radius 


R 


1-x 2 

friction-factor multiplier = 


s 


slip velocity ratio = (C /C T ) 


U 


rotor peripheral velocity 


w 


fluid velocity relative to rotor 


X 


quality E (m v /m T ) 


a 


void fraction = (A^/A^) 


B 


angle of fluid vector relative to tangent to 
rotor periphery 


B' 


angle of tangent to rotor-blade centerline relative 
to tangent to rotor periphery 
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A 



property difference 
deviation angle 
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e flow ratio = (Qb e /Q t p) 

p fluid density 

d) flow coefficient = (C /U) 

r m 

2 

ip work coefficient = (g Ah /U ) 

T co 

2 

ip’ head coefficient = (gAH Q /U ) 

SUBSCRIPTS 

0 total (static plus dynamic) or stagnation value 
of property 

1 plane at entrance to rotor 

2 plane at exit from rotor 

3 plane at entrance to pump 

be best-efficiency point 

f wall-friction component 

L liquid 

m meridional or radial component of velocity 

0 tangential component of velocity 

s flow-separation component 

sp single-phase component 

T total flow, liquid plus vapor 

tp two -phase flow 

th theoretical or ideal conditions 

v vapor 
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INTRODUCTION 



BACKGROUND 

In a nuclear power plant of the pressurized-water- 
reactor type (PWR) , sub-cooled water at high pressure and 
temperature is circulated through the reactor core by several 
main-coolant pumps, each in its own coolant leg. Safety 
analyses of postulated loss of coolant accidents (LOCA) 
require models of the flow through the main-coolant pumps 
for any one of several operating modes in both single- and 
two-phase flow. In an effort to more precisely predict the 
contribution of the main-coolant pumps to flow through the 
reactor core during a LOCA, there has been a continuing 
effort to improve these models, particularly for two-phase 
flow. 

A semi-empirical method for predicting the two-phase 
flow performance of a centrifugal pump was proposed by 
Mikielewicz and Wilson,^ and was applied to normal pump 
operation (forward rotation and flow) with promising results. 
Details of their procedure are summarized below. 

Starting with the Euler equation a theoretically- 
ideal head-versus-f low relationship is derived for 
both single- and two-phase flow. 

By examining the hydraulic losses in two-phase 
flow it is shown that for a given pump the relative 
pressure losses — the difference between theoretical 
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and actual pressure rise — are related to the relative 
pressure losses in single-phase flow by a function of 
void fraction, a property of the two-phase fluid. 

The ratio of relative pressure losses in two- 
phase flow to those in single-phase flow is defined 
as the "Head- loss ratio". 

The head-loss ratio is then calculated by 
combining the derived theoretical characteristics 
with experimental data for single- and two-phase flow. 

For first quadrant operation, good correlation is 
shown between the head-loss ratio and void fraction, 
and the result can be used to predict the two-phase 
characteristics of the pump. 

It is shown that the head-loss ratio for a given 
pump is related to other pumps by a function of pump 
geometry; further experimental work is required to 
determine this relationship. 

PURPOSE 

The purpose of this report is to determine if the 
procedure is applicable to the third quadrant of operation 
(reverse-flow and rotation). The steps of the procedure, as 
outlined, will be modified only as required by the different 
flow geometry, and the terminology used will be the same. 
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PROCEDURE 



DERIVATION OF THEORETICAL CHARACTERISTICS 

A centrifugal pump subjected to conditions such that its 
normal flow and rotation directions are reversed will operate 
in a manner similar to a radial-in-flow turbine. The Euler 
equation can be applied to determine the characteristics of 
the pump in this mode of operation: 

g c Ah Q = U lC01 - U 2 C 92 ^ 

where : 

Cg = fluid tangential velocity 

g c = constant in Newton's Law 

Ah Q = change in total enthalpy = h^-h^ 

U = rotor peripheral velocity 
subscript 1 = rotor entrance 
subscript 2 = rotor outlet 

The fluid tangential velocity at rotor outlet, C$ 2 ' 
could be significant because conditions of flow and rotor 
speed would be other than optimum and because, in general, 
a pump is not designed to be an efficient turbine. For 
these reasons the simplifying assumption of zero outlet 
swirl will not be made. This general form of the equation 
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is exact and will yield the actual enthalpy change across 
the pump, for any flow condition, if the proper values for 
the velocities are used. 

In the inlet region of the pump. Figure 1, precise 
flow angles and velocities are unknown and an assumption 
must be made relating C Q1 to the flow rate entering the 
scroll. Using this assumed value of the velocity tri- 

angle at the rotor entrance, as shown in Figure 2, can be 
solved for all velocities in terms of the flow rate. 

At the rotor outlet, flow geometry is specified by the 
design of what is normally the impeller inlet eye. A typical 
design with swept-back blading is shown in Figure 3. Reverse 
flow and rotation will produce a velocity diagram as shown in 
Figure 4, which can also be solved for all velocities in 
terms of the flow rate. 

If the reverse flow through the pump is idealized by 
making simplifying assumptions, the Euler equation can be used 
to approximate the theoretical head and flow characteristics 
of the pump. The following flow conditions are assumed. 

1. That the fluid tangential velocity at rotor entrance, 

Cg^, is equal in magnitude to the absolute fluid 
velocity, entering the scroll. 

2. That the relative flow angle at rotor outlet, & 2 > 
is constant for all flow rates and conditions. 

That the deviation angle at rotor outlet, & 2 ' is 

zero, that is, & 2 = & 2 * 

11 



3. 



PLANE 3 




FIGURE 1 - TYPICAL CENTRIFUGAL PUMP IN REVERSE-FLOW AND 
ROTATION 
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FIGURE 2 




U^ - rotor peripheral velocity 

- fluid absolute velocity 
Cg^- fluid tangential velocity 
C m ^~ fluid radial velocity 

- fluid relative velocity 
8^ - fluid relative flow angle 

- VELOCITY DIAGRAM AT IMPELLER INLET 
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FIGURE 3 - TYPICAL IMPELLER-EYE WITH REVERSE-FLOW AND 
ROTATION 
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FIGURE 4 



* 




- VELOCITY DIAGRAM AT IMPELLER OUTLET 
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With the above assumptions, equation 1 becomes 



Ah 



o th 



" u i c ei 



- U 2 C 02 



( 2 ) 



where 

subscript th s theoretical or ideal conditions 

✓ 

For the specific case of loss-less incompressible flow, 
the theoretical total enthalpy change is proportional to the 
change in theoretical total head.'*' 

9c 4t fe th = g ai tth 

change in total head = H ,-H^ 
gravitational acceleration 

and the expression for the theoretical total head character- 
istics of the pump is: 

9 AH oth = U 1 C 01 - °2 C 02 <3) 



where 

AH o = 

g = 
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Theoretical single-phase-f low characteristics 



For single-phase flow, the velocity triangles are solved 
by using the continuity equation and the assumptions made to 
idealize the flow. 

In the pump inlet region, it was assumed that 



where 

subscript 3 = plane at inlet to pump 



From continuity 



Q = 



C 3 A 3 



C ml A 



1 



where 

Q = volume flow per unit time 
A = area normal to the flow 

subscript m = meridional or radial component of velocity 



therefore , 



A 

^ 1 C mi a 



1 

3 



(4) 
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At the rotor outlet, from Figure 4, 



'92 



= U - JBi 

u 2 tan 8. 



Since 



U 



2 




r 



2 



r 



1 



where 



2 ... 
— = rotor-outlet-to-mlet radius-ratio 

r l 



and by continuity 




A 

C ml A 



1 

2 



then. 



C 



92 





ml 

tan 8. 



(5) 



By substituting equations 4 and 5, equation 3 becomes 



g AH othsp " U l C ml Aj " U 1 F^ [U 1 




C ml , 
tan 8 2 



where 



subscript sp = single-phase flow 
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Dividing by and defining the head coefficient, ij>^, as 



gAH 

n E T 

u: 



and the flow coefficient, <J>^, as 



_ C ml 



U. 



equation 3 becomes 



r 2 2 A 1 r 2 A 1 
♦ithsp = * 'a7 + A 2 tan B 2 1 h 



( 6 ) 



This expression is a linear function of ip ' versus <p , 
and is similar to the expression derived for first-quadrant 
operation . 



Theoretical two-phase- flow characteristics 

In two-phase separated flow it is assumed that the 
liquid and the vapor phases have separate mass flows and 
that their velocity components may also be different. In 
this case, equation 2 becomes 



U. 



^ithtp u 1 [ (1 “ x i )c 9Li +x i c e vi ] 



u: 



' [(1 " x 2 )C eL2 +x 2 C ev2 ] 
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where 



• • 

x = quality = 

and m = mass flow per unit time 
subscripts L = liquid 

T = total flow, liquid + vapor 
tp= two-phase flow 
v = vapor 



From the velocity triangles at rotor inlet and outlet 
and from the assumption that Cg^ = Cg, 



^lthtp 



1 A i ^1 

— [(l-x,)C T , t — + X-C , T-i] 
1 mLl Ag 1 mvl Ag 



U o C . O 

^2 C(1_X 2 ) (1 ~U 2 tan Bg 



) +x (i — siyJ )i 

; 2 u U„ tan B- ; J 



From continuity 






L 



m L m T (1-x) 

P L A L P L ( 1-a ) A T 



m 



Stlv 



v 



m T x 



P A 
v v 



p a A 
M v T 



where 

a = void fraction 
p = fluid density 
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and by defining the two-phase-flow coefficient, <J>^p, as 






Q, m_ 

_ *tp _ _T_ 



tp" A U p A u 



where 



P t p = mean density of two-phase flow 
= (l-a)p L +ap v 



Equation 2 becomes 






2 (1_Ct l ) P L1 



'1 a. 



vl 



+ X 2 ] - — {l- ^ .^P 2 [1 + 

X l ] tan V 1 



r ,, , 2 , (1 P L2 2, 1 

[(1-x.) + — x ]} 

2 a 2 p v2 2 



a _2 p v2 ■ 

(l-a 2 ) p l2 ' 



(7) 



By defining a two-phase-flow function, a. 



a = 



a 



'1-cr 




and a two-phase slip-velocity ratio, s. 
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The quality of the fluid can be expressed as 



as 

x = l+Ss 



By substituting these expressions, equation 7 becomes 



. , _ A 1 (l +a 1 ) ( 1+a i s i> _ r 2 2 

1 th tP a 3 tpl ( l+a^s 1 ) 2 r l 



[ 1 - 



4> 



t P 2 



(l+a 2 ) (l+a 2 s 2 ) 



tan Bo , , . ,2 

2 (l+a 2 s 2 ) 



-] 



(8) 



Applying the continuity equation to the two-phase flow 
coefficients at rotor inlet and outlet yields, 



p. , A, U. 

A _ A tpl 1 1 

V2 - *tpl P tp2 A 2 U 2 



which reduces equation 8 to 






- (— ) 2 [l-T^i^T!- f 



P t P l *1 “l. 



lthtp A^ tpl tpl tan 62 tp 2 p tp 2 A 2 U 2 



where 



( 1 +a) ( 1 +as ) 



tp ( 1 +as ) 2 



( 8 a) 
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Rearranging : 



. , - _/ 2.2 , ' 1 f , L 2 H tpl "1 tp2 , , 

^lthtp r^ A 3 tpl r;L P t p 2 a 2 tan B 2 ^tpl 

(9) 

If the flow becomes single-phase, equation 9 reduces to 
the single-phase case, equation 6, because f f 2 a nd 

P t p^/P t p2 all equal unity when the void fraction, a, 
equals zero. 

As derived, equation 9 is cumbersome and requires a 
knowledge of the relationships between inlet and outlet 
conditions to calculate values of the head coefficient for 
various conditions of void fraction and flow. By assuming 
that fluid properties at rotor inlet and outlet are equal, 
equation 9 reduces to 



^lthtp 




+ 



A i r2 A^ 

^Aj + A 2 tan 6 2 ] f tpl ^tpl 



( 10 ) 



The assumption is valid in this case because the moment of 

momentum at rotor outlet, U C _, is very small in comparison 

to U.C.., and the outlet properties appear in the term derived 
1 01 

from U 2 C 0 2 ; t ^ lus the e ff ect °f changes in fluid properties 
from inlet to outlet is minimal. 
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The two-phase flow function, f 



was defined by 



tp' 



equation 8a as 



tp 



(1+a) (1+as ) 
(1+as) 2 



where 



a = 



(— ) — 



s = 



v 



For a steam-water mixture, liquid and vapor densities 

are given by steam tables for various pressures. The slip- 

velocity ratio, s, for two-phase flow in a horizontal pipe, 

2 

was correlated as a function of pressure and void fraction 

as shown in Figure 5. A correlation more applicable to two- 

phase flow in centrifugal pumps was not available; therefore, 

the relationship described, though for vastly different 

conditions of flow, was applied. 

Values of the two-phase flow function, f , can be 

calculated as a function of void fraction, a, for various 

pressures, as shown in Figure 6. 

The theoretical characteristics for third-quadrant 

operation of a centrifugal pump in both single- and two- 

phase flow can be represented as shown in Figure 7. In 

single-phase flow, f is unity; in two-phase flow, f fc ^ is 

determined by pressure and void fraction. 
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FIGURE 5 - SLIP-VELOCITY RATIO, s, VERSUS PRESSURE 



(from reference 2) 
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1.201 



1. 18- 




VOID FRACTION, a 



FIGURE 6 - TWO-PHASE FLOW FUNCTION, f , VERSUS VOID 

tp' 

FRACTION, a, FOR 200 AND 500 psia 
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n t h 





f 



t P i 



<p 



i 



FIGURE 7 - THEORETICAL HEAD-VERSUS-FLOW CHARACTERISTICS 
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FLOW LOSSES 



At a given rate of flow through a centrifugal pump, 

the difference between the measured head and the theoretically 

ideal head is the result of flow losses generated within the 

pump, and changes in the exit deviation angle , 6 2 , from the 

assumed ideal direction. For third-quadrant operation, in 

single-phase flow, this difference is shown in Figure 8. 

The curve of actual head versus flow is typical for this 
3 4 

quadrant, ' and the theoretical relationship is approxi- 
mated from the previously derived expression (equation 6) . 

The actual head, AH , is a function of the flow-rate, Q, 

o 

squared;^ therefore, with flow increasing from zero, the 
losses decrease as the square of the flow, pass through a 
minimum, then increase as the square of the flow. The flow- 
rate corresponding to the point of minimum losses is de- 
fined as the best efficiency flow, , for third-quadrant 
operation . 

In the first-quadrant of operation, the flow losses 
a 

behave in^similar manner and the expressions derived to 
relate the losses in two-phase flow to the losses in single- 
phase flow may be applied to the third-quadrant. The 
constants and flow-rate ratio are those applicable to the 
third quadrant. 
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FIGURE 8 - TYPICAL THEORETICAL AND ACTUAL PUMP 
CHARACTERISTICS 
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The flow losses are expressed as the sum of wall-friction 



losses and flow-separation losses, 

Ap = Ap j. + Ap 
^o ^of os 



( 11 ) 



where 

Ap Q = loss of total, or stagnation, pressure 
subscripts f = friction component 

s = flow-separation component 



The wall-friction component of the losses in two-phase 
flow is given by 



Ap . , = R K 

^otpf 



^P)2 



Q 



tp 



( 12 ) 



where 

5 1-x 2 

R = friction factor multiplier [j- 3 ^] 

t 

K = constant at high Reynolds numbers; includes pump 
geometry, the single-phase friction factor, and 
liquid density 

By approximating the flow-separation losses as proportional 

to the square of the difference between the existing flow 

6 

and the flow at best efficiency point, Q^/ the two-phase 
flow-separation losses become 
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(13) 



Ap 



otps 



= K. 



•tp 



[ (1-ct) ( 



1+a 

1+sa 



) - e] 



where 



= constant, depends on flow regime 
s 5 Q be /Q tp 



Adding these components yields the expression for the 
total losses in two-phase flow, Ap Qt p, 



4p otp ■ R K <pf> 2 Q tp + K 4 Otp'* 1 -”) ' e ’ 2 < 14 > 



substituting 



p tp _ r l—oi -I f l+a , 

P L L l-x J 1+sa 

The total losses in two-phase flow become 



Ap 



otp 



- V<in «> 2 + 



k 6 1(1 '“> ( irfc' ) ' El2} Q t P 



(15) 



where 

Kg and Kg are constants which include geometric and 
flow-regime factors 

for single-phase flow, equation 15 reduces to 
Ap osp - V 1 + K 6 (1 - e)2 l °L 
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The comparison between single-phase and two-phase flow 
is made at like values of flow-rate; therefore, the ratio 
of two-phase to single-phase losses is 



A P, 

Ap 



,1+a >2 . r , i \ ,1+a > -.2 

ot£ = ‘lTi? 1 + K 6 1(1 - a) ( ITS¥ ) ~ 



osp 



1 + K (1-e) 

D 



(17) 



HEAD-LOSS RATIO 

The total losses, Ap , can be expressed as the difference 

o 

between the actual and the theoretical heads at a given flow- 
rate , 




(AH 

o 



AH 
o th 




(18) 



In terms of head, equation 17 becomes 



AH 



AH 



o tp 



o sp 



AH 



o thtp _ 



AH 



o thsp 




( l±a_) 

l l+sa ; 



2tK 6 [(1 - a)( irfi- ) - £l2 

1+Kg (1-e) 2 



(19) 



substituting 



tp 



, 1— x ^ ,1+a 
\ 1 ^ ) \ * 



'1-a l+sa ; 



and x 



sa 

1+sa 
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Equation 19 becomes 



AH 



f 1+a 



,1+a 



otp AH othtp _ ( l+sa ) + K 6 [(1 a)( l+sa ) ~ e] 



AH_ - AH .. 
osp othsp 



(1-a) (1+a) [1+K g (1-e ) z ] 



( 20 ) 



For the first quadrant, the head-loss ratio, H* , was 
defined as 



H* 



AH 



othtp 



AH 



othsp 



AH 

AH 



otp 

osp 



To prevent the introduction of an unnecessary change, 
the first quadrant definition will be used in the third 
quadrant as it is numerically equal to the ratio in 
equation 20. Thus, for third quadrant operation 



H * ° - AH° tp = £[ <£§I>' V el 



( 21 ) 



othsp 



osp 



In this functional relationship: (1+a/l+sa) is a 

function of void fraction and pressure; Kg includes the 
geometry of a given pump and the flow regime of the two- 
phase fluid (a function of void fraction) ; and e is the 

i 

ratio of best efficiency flow in third-quadrant operation 
to the existing flow. For a given pump, operating at a 
particular pressure, the ratio of two-phase to single-phase 
losses at like values of e is therefore principally a 
function of the void fraction of the two-phase fluid. 
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Substituting the definition for the head coefficient 









gAH 



U 



the head-loss ratio can be expressed as 



H * = ^ ht P 

\p! u - \p ' 

r thsp r sp 



( 22 ) 



CALCULATION OF THE HEAD-LOSS RATIO 

The head-loss ratio as a function of void fraction for 
third-quadrant operation of a particular centrifugal pump 
can be constructed using the derived theoretical character- 
istics and the actual characteristics obtained from experi- 
mental data. The relationship between H* and a for this 
pump can then be used to predict the third-quadrant, two- 
phase performance of other pumps as single-phase character- 
istics are usually known or can be easily determined. 

D.J. Olson conducted experiments on a semi-scale pump 
of known geometry, ^ in single- and two-phase flow, and in 
the operating quadrants of interest. Using the results of 
his experiments in the third quadrant, the actual character- 
istics of the pump can be calculated. Information obtained 

g 

from the pump manufacturer provided the required constants 
for calculation of the theoretical characteristics for the 
pump using equation 10. 
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The functional relationship between head-loss ratio 
and void fraction for this particular pump is developed as 
described in the following sections. Details of Olson's 
apparatus and test conditions are shown in Appendix D. 

Theoretical head-coefficients 

Details of the semi-scale pump used by Olson are shown 
in Figures 9 and 10. With proper values substituted for the 
constants, equation 10 becomes 

^ th = -0.1896 + 9.155 f tp] _ <i> 1 

As previously described, this equation may be used to 
determine the theoretical characteristics in both single- 
and two-phase flow. 

Actual single-phase head-coefficient 

A non-dimensional head-flow curve for the semi-scale 
pump can be constructed, as shown in Figure 11, using the 
single-phase data. Curve-fitting the points by regression, 
using the method of least squares, this curve is accurately 
approximated by the expression 

4>| sp =101.28 ^ 1 1 * 8756 + 0.4 
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FIGURE 9 




- DETAILS OF THE SEMI -SCALE PUMP 
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FIGURE 10 - DETAILS OF THE SEMI-SCALE PUMP IMPELLER EYE 
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0.05 0.1 0.15 0.2 




FIGURE 11 - HEAD- VERSUS -FLOW CHARACTERISTICS OF SEMI-SCALE 
PUMP IN SINGLE-PHASE, THIRD-QUADRANT OPERATION 
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Head-loss ratio 



The head-loss ratio is calculated for each test 
conducted in two-phase flow, third-quadrant operation. 

RESULTS 

The results of these calculations are shown in Table 1, 
and the relationship between head-loss ratio and void 
fraction for the semi-scale pump is shown in Figure 12. 
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a 



1 



FIGURE 12 - HEAD-LOSS RATIO, H* , VERSUS INLET VOID FRACTION, 



40 



TABLE 1 



CALCULATION OF HEAD-LOSS RATIO, H* 



INLET 

VOID 

FRACTION 


f tpl 


♦tpi 


Hthtp 


*ltp 


^lthsp 


^lsp 


H* 


a l 
















0.461 


1.022 


0.030 


0.095 


0.919 


0.089 


0.545 


1 . 


0.464 


1.015 


0.018 


-0.022 


0.876 


-0.025 


0.454 


1 . 


0.466 


1.022 


0.010 


-0.094 


0.816 


-0.096 


0.419 


1 . 


0 


1.0 


0.009 


-0.102 


0.469 


-0.102 


0.417 


1 . 


0 


1.0 


0.010 


-0.101 


0.481 


-0.101 


0.417 


1 . 


0 


1.0 


0.024 


0.027 


0.570 


0.027 


0.491 


1 . 


0.241 


1.015 


0.044 


0.215 


0.824 


0.209 


0.683 


1 . 


0.116 


1.009 


0.047 


0.247 


0.811 


0.244 


0.732 


1 . 


0.245 


1.007 


0.044 


0.214 


0.848 


0.211 


0.686 


1 . 


0.056 


1.004 


0.025 


0.043 


0.629 


0.042 


0.502 


1 . 


0.512 


1.067 


0.037 


0.168 


1.004 


0.146 


0.605 


1 . 


0.537 


1.023 


0.022 


0.012 


0.956 


0.008 


0.476 


2. 


0.631 


1.033 


0.015 


-0.048 


1.184 


-0.052 


0.438 


2 . 


0.401 


1.014 


0.016 


-0.043 


0.788 


-0.046 


0.442 


1 . 



81 

88 

77 

10 

12 

17 

28 

16 

33 

27 

82 

02 

51 

70 



41 



DISCUSSION 



Figure 12 shows that within the range of variables of 
Olson's experiments, the head- loss ratio of the semi-scale 
pump in third-quadrant operation correlates very well as a 
function of inlet void fraction. An average curve drawn 
through the points can be used to predict the two-phase 
performance at other values of void fraction, flow, and 
pressure, up to the maximum void fraction tested. 

At low values of void fraction (up to 20%) the flow is 
largely homogeneous and the losses in two-phase flow are 
approximately equal to those in single-phase. As void 
fraction increases above 20%, the flow becomes separated 
and the relative losses increase sharply, until at the 
maximum void fraction tested (63%) , they are more than 
twice the single-phase losses. By the definition of the 
head-loss ratio, its value is unity at a void fraction of 
100%; the characteristics between 63tand 100% void-fraction 
are unknown, but it is reasonable to assume that the head- 
loss ratio will continue to increase then drop sharply to 
unity at the very high values of void fraction. 

The effect of pressure on the correlation of the head- 
loss ratio is not apparent for this limited amount of data. 
Two series of two-phase tests were conducted, one at 200 psia 
system pressure, the other at 500 psia. The functional 
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relationship for the head-loss ratio, equation 21, indicated 
that H* was dependent on void fraction, a, and pressure 
(for a given pump) . in this equation, system pressure 
affects the two-phase flow function, a, through its effect 
on density and the slip-velocity ratio, s. Both of these 
effects are overshadowed by changes in void fraction, a. 

OTHER RESULTS 

The head-loss ratio was calculated using expressions 
other than equation 10 for the theoretical two-phase head 
coefficient. 

The experimental data for the fluid properties at 
pump-outlet were used to evaluate equation 9 and determine 
the head- loss ratio. Results are shown in Table El and 
Figure El. These results are not significantly different 
from the results presented in Figure 12, indicating that 
the effect of the outlet tangential velocity, Cq2’ 
minimal, as assumed. Equation 10 can be applied to any 
pump, without prior knowledge of the relationships between 
fluid properties at inlet and exit; therefore, it is much 
more useful as a means of predicting two-phase characteristics. 

Results obtained by assuming zero outlet-swirl are 
shown in Figure E2, and also show good correlation between 
H* and a, but indicate much higher losses. The procedure is 
simplified by this assumption, and further experiments may 
prove it to be viable for certain pump designs, particularly 
those of low specific speed. 
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APPLICATION OF THE PROCEDURE 



The two-phase, third-quadrant head-ver sus-f low 
characteristics of a centrifugal pump can be predicted, for 
various values of void fraction and pressure, using the pro- 
cedure outlined below. The following information concerning 
the pump must be known . 

Flow area at the pump discharge flange, A^ • 

Flow area of the outer radius of the impeller, A^. 

Flow area at the impeller eye, A^. 

The impeller radius ratio, r^/r^. 

Blade angle at impeller eye, B^, as defined by 
Figure 3. 

Sufficient single-phase data to construct a plot 
of head-coefficient versus flow-coefficient. 

The actual single-phase, and the theoretical single- 
and two-phase head-versus-f low characteristics are determined 
by constructing a graph as shown in Figure 13. The three 

head-coefficients, ^ithsp' an< ^ ^lthtp' Can then ke 

found for particular values of flow-coefficient, <J>^, and 
two-phase flow function, f^ 

At the desired value of void fraction, a^, and system 
pressure, the head- loss ratio, H*, is determined from 
Figure 12, and the two-phase flow function, f calculated. 
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FIGURE 13 - ACTUAL SINGLE-PHASE AND THEORETICAL SINGLE - 
AND TWO-PHASE PERFORMANCE 
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The predicted values of actual two-phase head coeffi- 
cient, , are calculated as a function of flow-coefficient, 

4>^, from 

*itp ‘ *itht P - H * [ *iths P -*i S p ] 

Experiments on pumps of different geometries will be 
required to determine the effect of geometrical factors such 
as specific speed, N , on the functional relationship of the 
head- loss ratio, H* (equation 21). Until such data becomes 
available, this procedure may be used to determine approximate 
two-phase characteristics. 
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APPENDIX A 



DETAILS OF THE DERIVATION OF THE THEORETICAL 
TWO-PHASE HEAD COEFFICIENT 



In two-phase separated flow, the theoretical two-phase 
head coefficient is given by 



^lthtp U.^ (1 " X l )C 6Ll +X l C evl ]_ u 2 [ (1 X 2 )C 0L2 +X 2 C 0v2 ] 



(Al) 



In general, velocities of the two phases are given by 



m T 



m T (l-x) 



CL ’ P L A L " C L (1 - a,A T 



m 



v 



m T x 



C = s 7\ 

v p A pa A m 
K v v v T 



Assuming Cg 1 = C^ 



m T (1_X l ) 

C 0 Ll = C L3 = p L1 d-a 1 ) a 3 



= C 



m T X 1 



'0vl v3 P v ^ a^ A^ 



By continuity 



m T = constant 



SmLl P Ll (1 ~ a l ) A 1 

1 - x. 



C , p . a, A. 
mvl vl 1 1 

x. 
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Then 



_ C mLl p Ll (1 ~ a l )A l (1 “ x l ) 



0L1 



(l-x 1 ) 



p Ll (1 '“l ,A 3 



and 



0L1 



C — 
mLl A 3 



Similarly 



C = c — — 

0vl mvl A 3 



Substituting into equation A1 



C T , A n C 1 A. 

I/, > = (1— x ) rnkl __1 i. 

dthtp u 1 ; a 3 1 u 3 a 3 



(i * 2 )U 2 C 0L2 x 2 U 2 C 0 v2 i 
l o + O J 



U. 



U, 



From the general expressions for velocity 
m T (1-x) 

C mLl “ p Ll ( 1-a^) A^ 

and 



C mvl p vl a 1 



(A2) 
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Substituting these expressions, the first two terms 



on the right hand side of equation A2 become 

2 • 2 ♦ 

(l-x 1 ) m T A 1 x x m T 

M * A- p a n UtA. A- 

U 1 P L1 ( 1-a i) A i 3 vl 1113 

Rearranging and multiplying by 



(l-a 1 ) 

(l-a 1 ) 



P L1 + °1 P vl 



P L1 + “l P vl 



where 



(1-a) p T + ap = p . 

' L v tp 



They become 



m A (1-x, ) 2 x 2 

p tpl *1 *3 ‘ p Ll U-a7T + P L1 + “1 P vl )! 

Substituting the definition of the two-phase flow coefficient 



m„ 



<P 



tpv P tpl A 1 U 1 



the two terms become 



a 7 ♦tpl [(1 " x l } + p^JO^iJ 



2 (! x x ) ct 1 P vl ^ x i P L1 (1 ~ ct i ) + x 2j 



p vl a l 
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by definition, the two-phase flow function, a, is 



a = 



r— i __z 
W p L 



and the quality, x, becomes 



x = 



as 

1+as 



where 



s = 



substituting these expressions, the first two terms become 



A, 2a i s i (a,s ) 

*tol tl+a l ] [1 -TTi-i- + — ^ 

A 3 tpl 1 1 a l 1 l+2a 1 s 1 +(a 1 s 1 ) 



2 < 1+ l> ' 



performing the indicated operations and rearranging yields 



A i , 

AJ *tpl 



(l+a 1 ) (l+a 1 s 1 ) 
(l+a 1 s 1 ) 2 



Substituting this expression, equation A2 becomes 



r. 



= IT- 4>_i f, 



lthtp A^ tpl tpl 



r (l' X 2 )U 2 C 0L2 1 x 2 U 2 C 0v2 1 

— L o " i * o J 



U, 



U 
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where 



f 



tpl 



(l+a ] _) (1+a.^ 2 ) 
(l+a 1 s 1 ) 2 



From the outlet velocity triangle. Figure 4 , the last two 
terms in equation A2 reduce to 



U 



f tp2 *t P 2 , , r _2. 2 

tan Q 2 1 v r ; 



using the same procedure applied to the first two terms. 
Equation A2 becomes 



^lthtp A„ ^tpl ^tpl 



2 2 

-(zA) [i 



f tp2 ^ tp2 , 
tan 6 2 1 



The two-phase flow coefficients, , at inlet and outlet 
can be related by the continuity equation. 



^ t P 2 



^tpl 



p tpl 

p t P 2 





Substituting this expression and rearranging yields 



♦ithtp - -<^> 2 + & £ t Pl 



+ h Pt p . 1 tk ffc p 2 - ] ^ 

r x p tp2 A 2 tan 6 2 tpl 



(A3) 
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APPENDIX B 



DETAILS OF THE CALCULATION OF HEAD-LOSS RATIO, H* 

The head-loss ratio, H* , for the various values of void 
fraction, a, was calculated using the relationships previously 
presented. The data from the experiments conducted by Olson 
are included in Appendix C. Values for fluid densities 
were interpolated from reference 9. 

Table B1 lists the single-phase characteristics used 
to determine the single-phase head coefficient, , as a 

function of flow coefficient, Tables B2 through B5 

present the calculation of the head-loss ratio, H* . 

The test number column refers to the experimental data. 
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TABLE B1 



CALCULATION OF THE SINGLE-PHASE HEAD COEFFICIENT, 





AS A FUNCTION OF FLOW 


COEFFICIENT 


1 / ^ ^ 




SINGLE- 

PHASE 

TEST 

NUMBER 


gAH 

2 2 
(ft /sec) 


U 1 

(ft/sec) 


C ml 

(ft/sec) 


lsp 7 
(gAH/u 1 ) 


*1 

(C ml /U 


78 


1330 


13.2 


3.188 


7.64 


0.24 


82 


1459 


14.2 


3.275 


7.22 


0.23 


86 


164.22 


13.4 


0.75 


0.91 


0.06 


87 


570 


13.5 


1.962 


3.15 


0.15 


88 


763 


13.5 


2.308 


4.22 


0.17 


89 


1346 


13.6 


3.188 


7.40 


0.23 


92 


74 


13.6 


0 


0.4 


0 


83 


4176 


20.8 


3.635 


9.63 


0.17 


75 


605 


27.1 


1.255 


0.82 


0.05 


76 


960 


27.2 


2.121 


1.30 


0.08 


77 


1513 


27.1 


3.231 


2.06 


0.12 


91 


299.5 


27.04 


0 


0.41 


0 


95 


2380 


53.78 


2.568 


0.82 


0.05 


98 


1204 


54.15 


0 


0.41 


0 


150 


4405 


54.25 


4.95 


1.50 


0.09 


106 


3130 


74.6 


1.414 


0.56 


0.02 


107 


3916 


74.7 


2.640 


0.70 


0.04 


108 


4634 


74.6 


3.534 


0.83 


0.05 


149 


7065 


74.5 


5.756 


1.27 


0.08 


147 


7467 


90.9 


5.02 


0.90 


0.055 


148 


9792 


91.2 


6.492 


1.18 


0.07 


152 


9982 


92.0 


6.65 


1.18 


0.07 
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TABLE B2 



CALCULATION OF TWO-PHASE FLOW FUNCTION, f , 
-LE. 



AND TWO-PHASE DENSITY, P , AT PUMP INLET 



TWO- 


VAPOR 


LIQUID 


VOID 








PHASE 

TEST 

NUMBER 


DENSITY 

Pvl 3 
lb m /ft 


DENSITY 

P]L1 3 
lb n/ ft 


FRACTION 

a l 


f t P l 


p tpl 

lb /ft 
m 


149 


0.4377 


54.38 


0.548 




1.0299 


24.82 


150 


0.4443 


54.33 


0.461 




1.0218 


29.49 


151 


0.4508 


54.28 


0.464 




1.0151 


29.30 


152 


0.4354 


54.40 


0.466 




1.0218 


29.25 


153 


0.4354 


54.40 


0** 




1.0 


54.40 


154 


0.4420 


54.35 


0** 




1.0 


54.35 


155 


0.4309 


54.43 


0** 




1.0 


54.43 


156 

* 


0.4288 


54.45 


0.241 




1.0150 


41.43 


157 


1.1002 


50.53 


0.116 




1.0086 


44.80 


158 


1.1066 


50.51 


0.245 




1.0069 


38.41 


159 


1.1042 


50.52 


0.056 




1.0039 


47.75 


160 


1.0932 


50.57 


0.512 




1.0671 


25.24 


161 


1.1156 


50.46 


0.537 




1.0232 


23.96 


162 


1.1156 


50.46 


0.631 




1.0327 


19.32 


163 


1.1453 


50.33 


0.401 




1.0143 


30.61 


164 


1.1202 


50.44 


0.647 




1.0348 


18.53 


*Tests 


149 through 


156 were 


conducted 


. at 


nominal 


system 


pressure of 200 psia, from : 


Figure 5, 


slip velocity ratio, 


s = 2.8 


. Tests 157 


through 


164 were 


conducted at 


nominal 


system 


pressure of 


500 psia 


, s = 2 . 









* * 



Negative void fractions in data were taken to be zero 
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TABLE B3 



TWO- 

PHASE 

TEST 

NUMBER 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 



CALCULATION OF TWO-PHASE FLOW FUNCTION, f , 



AND TWO-PHASE DENSITY, p , AT PUMP OUTLET 



VAPOR LIQUID VOID 



DENSITY 

Pv2 3 
lb n/ ft 


DENSITY 

P L2 

lb r/ ft 


FRACTION 

a 2 


f tp2 


p tp2 

lb /ft 
m 


0.4261 


54.47 


0.840 


1.1070 


9.07 


0.4113 


54.60 


0.867 


1.1234 


7.62 


0.4197 


54.52 


0.781 


1.0767 


12.27 


0.4197 


54.52 


0.558 


1.0298 


24.33 


0.4009 


54.68 


0.008 


1.0002 


54.25 


0.3967 


54.71 


0.007 


1.0002 


54.33 


0.3946 


54.73 


0.730 


1.0568 


15.07 


0.3864 


54.80 


0.789 


1.0742 


11.87 


1.0514 


50.76 


0.702 


1.0405 


15.86 


1.0646 


50.70 


0.756 


1.0509 


13.18 


1.0624 


50.71 


0.428 


1.0147 


29.46 


1.0646 


50.70 


0.696 


1.0400 


16.15 


1.0844 


50.61 


0.912 


1.1065 


5.44 


1.0844 


50.61 


0.988 


1.0860 


1.68 


1.1173 


50.46 


0.744 


1.0505 


13.75 


1.0998 


50.54 


0.834 


1.0736 


9.31 



(1+a) ( 1+as 2 ) 
(1+as) 2 




(1-a) P L + ap v 
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TABLE B4 



CALCULATION OF 


TWO- PHASE 


FLOW COEFFICIENT, <|> tpl 


t 




AND TWO- 


PHASE HEAD 


COEFFICIENT, 


*itp 




TWO- 

PHASE 

TEST 

NUMBER 


FLOW 
RATE, Q 

(gal/min) 


ROTOR 
SPEED, N 

(rev/min) 


^tpl 

(0.951 Q/N) 


HEAD 

AH 

o 

(ft) 


Ht P 

(28,200 

AH^ 

°) 

2 ' 

N 


149 


<J\ 

• 

in 

00 

I 


0 


— 


40.2 


— 


150 


-51.5 


-1608 


0.030 


84.3 


0.919 


151 


-30.3 


-1600 


0.018 


79.6 


0.876 


152 


-17.3 


-1607 


0.010 


74.8 


0.816 


153 


-16.2 


-1604 


0.009 


42.8 


0.469 


154 


-16.4 


-1602 


0.010 


43.8 


0.481 


155 


-39.8 


-1600 


0.024 


51.8 


0.570 


156 


-73.4 


-1603 


0.044 


75.1 


0.824 


157 


-80.0 


-1607 


0.047 


74.3 


0.811 


158 


-73.8 


-1603 


0.044 


77.3 


0.848 


159 


-42.6 


-1600 


0.025 


57.1 


0.629 


160 


-61.4 


-1594 


0.037 


90.5 


1.004 


161 


-36.2 


-1597 


0.022 


86.5 


0.956 


162 


-25.3 


-1600 


0.015 


107.5 


1.184 


163 


-26.5 


-1603 


0.016 


71.8 


0.788 


164 


-112.9 


0 


— 


76.4 


— 
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TABLE B5 



CALCULATION OF HEAD-LOSS RATIO, H* 



TWO- 

PHASE 

TEST 

NUMBER 


^lthtp 


♦it P 


^lthsp 


4) ' 

lsp 


H* 


149 


— 


— 


— 


— 


— 


150 


0.095 


0.919 


0.089 


0.545 


1.81 


151 


-0.022 


0.876 


-0.025 


0.454 


1.88 


152 


-0.094 


0.816 


-0.096 


0.419 


1.77 


153 


-0.102 


0.469 


-0.102 


0.417 


1.10 


154 


-0.101 


0.481 


-0.101 


0.417 


1.12 


155 


0.027 


0.570 


0.027 


0.491 


1.17 


156 


0.215 


C.824 


0.209 


0.683 


1.28 


157 


0.247 


0.811 


0.244 


0.732 


1.16 


158 


0.214 


0.848 


0.211 


0.686 


1.33 


159 


0.043 


0.629 


0.042 


0.502 


1.27 


160 


0.168 


1.004 


0.146 


0.605 


1.82 


161 


0.012 


0.956 


0.008 


0.476 


2.02 


162 


-0.048 


1.184 


-0.052 


0.438 


2.51 


163 

164 


-0.043 


0.788 


-0.046 


0.442 


1.70 




n t h = 


-0.1896 


+ 9.155 


f tpi n 






*isp = 


101.28 <j> 


1.8756 

1 


0.4 





H * = ^ithtp ^itp 



INLET 

VOID 

FRACTION 



0.461 

0.464 

0.466 

0 

0 

0 

0.241 

0.116 

0.245 

0.056 

0.512 

0.537 

0.631 

0.401 
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APPENDIX C 



EXPERIMENTAL DATA 

Tables Cl and C2 are reproduced from reference 7. 

Explanatory notes for single-phase tests: 

a. Based on turbine flowmeter measurements. 

b. Positive pump AP represents higher pressure 
at discharge nozzle. 

c. Based on fluid density measurement at pump 
inlet. 

d. Fluid density at pump inlet (suction for 
normal operation, discharge for reverse flow) . 

e. Inlet and outlet refer to locations with 
respect to direction of flow. 

f. Flow rate obtained from venturi flowmeter 
in test loop. 

Explanatory notes for two-phase tests: 

a. Inlet flow based on loop mass flow rate and 
test section inlet (suction or discharge) 
fluid density. 

b. Positive AP indicates higher pressure at 
pump discharge. 

c. Based on inlet conditions. 

d. Inlet and outlet refer to the direction of 
flow. 
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e. Calculated using measured density and 
assumed saturation conditions. 

f. Measured directly. 
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RESULTS OBTAINED FROM TEST SECTION MEASUREMENTS DURING SINGLE-PHASE PLOW PUMP TESTS 
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APPENDIX D 



DESCRIPTION OF EXPERIMENTS 

The semi-scale pump used in the experiments had the 
following characteristics: 

Rated flow - 180 gpm 
Rated head - 192 ft. 

Specific speed - 926 
Rated torque - 417.6 in-lb 
Details of the pump are shown in Figures 9 and 10. 

The test-loop used for single-phase tests is shown in 
Figure Dl. The two-phase test-loop is shown in Figure D2 . 
Superheated steam and/or saturated water were supplied by 
an oil-fired boiler. 

Flow conditions imposed for third quadrant tests are 
as follows. 

Single-phase tests 

Impeller speed - up to minus 2700 rpm 
Flow rate - 0 to minus 207 gpm 
Pressure at pump inlet - 60.7 to 216.6 psig 
Pump head - 5.1 to 310.0 ft of fluid 
Two-phase tests 

Impeller speed - 0 and minus 1600 rpm 

Flow rate - minus 16.2 to minus 112.9 gpm 

System pressure - 200 and 500 psia 

Pump head - 40.2 to 107.5 ft of fluid 
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P Pressure Top 

Differential Pressure Transducer 
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FIGURE Dl - SINGLE-PHASE TEST LOOP 
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(from reference 7) 




FIGURE D2 - TWO-PHASE TEST LOOP (from reference 7) 
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APPENDIX E 



OTHER RESULTS 

This section summarizes, for purposes of comparison, 
the correlation of H* vs a obtained using alternate methods 
to compute the theoretical characteristics. 

The experimental data were used to evaluate equation 9, 
thereby making it semi-empirical rather than theoretical. 

The head coefficients obtained in this manner were used to 
calculate the head-loss ratio. Results are shown in Table El 
and Figure El. 



ZERO SWIRL AT OUTLET 

The theoretical characteristics can be derived assuming 
that outlet swirl is negligible. In this case the expression 
for the theoretical head coefficient, is 



^lth A 3 f tpl *1 



Results are shown in Table E2, and the correlation of 
H* versus ^ is shown in Figure E2. 
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FIGURE El - HEAD-LOSS RATIO VERSUS INLET VOID-FRACTION 
USING EQUATION 9 FOR 'i'{ thsp 
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TABLE El 



CALCULATION OF THE HEAD-LOSS RATIO , H* , USING 
EXPERIMENTAL DATA TO EVALUATE EQUATION 9 



TWO- 

PHASE 

TEST 

NUMBER 


( ^lthtp ) se 


H* 


a l 


150 


0.199 


1.581 


0.461 


151 


0.006 


1.816 


0.464 


152 


-0.092 


1.763 


0.466 


153 


-0.102 


1.100 


0 


154 


-0.101 


1.124 


0 


155 


0.096 


1.023 


0 


156 


0.336 


1.030 


0.241 


157 


0.341 


0.964 


0.116 


158 


0.306 


1.142 


0.245 


159 


0.059 


1.238 


0.056 


160 


0.189 


1.776 


0.512 


161 


0.097 


1.835 


0.537 


162 


0.128 


2.156 


0.631 


163 


-0.022 


1.660 


0.401 


^lthtp 


= -0.1896 + (8.1342 


f t . + 1.021 
tpl 


P tpl f 
°t P 2 t ? 2 



Values for f f tp2 ' p tpl' and p tp2 were obtained from 

Appendix B. 



t 



where subscript se 



semi-empirical 
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FIGURE E2 - HEAD-LOSS RATIO VERSUS INLET VOID FRACTION FOR 
ZERO-SWIRL AT OUTLET, C 02 = 0. 



70 



TABLE E2 





CALCULATION OF THE 


HEAD-LOSS RATIO, 


H*, 






ASSUMING ZERO SWIRL 


AT ROTOR OUTLET, 




= 0 


TWO- 

PHASE 

TEST 

NUMBER 


^Ispth 


'Htpth 




H* 


150 


0.248 


0.253 




2.24 


151 


0.146 


0.148 




2.36 


152 


0.083 


0.085 




2.18 


153 


0.078 


0.078 




1.15 


154 


0.079 


0.079 




1.19 


155 


0.193 


0.193 




1.27 


156 


0.354 


0.359 




1.41 


157 


0.385 


0.388 




1.22 


158 


0.356 


0.358 




1.48 


159 


0.206 


0.207 




1.43 


160 


0.298 


0.318 




2.23 


161 


0.175 


0.179 




2.58 


162 


0.122 


0.126 




3.35 


163 


0.128 


0.130 




2.10 
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